m 



^ 



> 

1^ 



New Regions in the NMSSM with a 125 GeV Higgs 



Marcin Badziakj^ Marek OlechowskQ Stefan Pokorsk|^ 

Institute of Theoretical Physics, Faculty of Physics, University of Warsaw 
ul. Hoza 69, PL-00-681 Warsaw, Poland 



O 
(N 

^ Abstract 

< 

Q\ It is pointed out that mixing effects in the CP-even scalar sector of the NMSSM can give 

^~~* 6-8 GeV correction to the SM-hke Higgs mass in moderate or large tan/3 regions with a small 

""I value of the singlet-higgs-higgs superfields coupling A ~ (9(0.1). This effect comes mainly from 

P^ the mixing of the SM-like Higgs with lighter singlet. In the same parameter range, the mixing 

o I of the heavy doublet Higgs with the singlet may strongly modify the couplings of the singlet- 



^ like and the 125 GeV scalars. Firstly, the LEP bounds on a light singlet can be evaded for 
a large range of its masses. Secondly, the decay rates of both scalars can show a variety of 
interesting patterns, depending on the lightest scalar mass. In particular, a striking signature 
of this mechanism can be a light scalar with strongly suppressed (enhanced) branching ratios 
cn to bb {gg, cc, '-ff) as compared to the SM Higgs with the same mass. The 77 decay channel 
i/-^ is particularly promising for the search of such a scalar at the LHC. The 125 GeV scalar can, 
■^ thus, be accommodated with substantially smaller than in the MSSM radiative corrections from 
5 the stop loops (and consequently, with lighter stops) also for moderate or large tan/3, with the 
^-H mixing effects replacing the standard NMSSM mechanism of increasing the tree level Higgs 
>> mass in the low tan/3 and large A regime, and with clear experimental signatures of such a 
^ mechanism. 
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1 Introduction 

The discovery of a SM-like Higgs particle has recently been announced by the LHC experiments 
[H, |2]. Although its properties such as the production and decay rates into different channels 
still remain very uncertain [3[ H], its mass is established to be around 125 GeV, with only a 
couple of GeV uncertainty, and this puts new constraints on the BSM models. In the minimal 
supersymmetric SM (MSSM) the Higgs tree-level quartic coupling is given by the electroweak 
gauge coupling, so that the theory predicts a tree-level upper bound for the Higgs mass to 
be equal Mz- It is well known that loop corrections, mainly from the top-stop loop, can 
significantly raise the Higgs mass in the MSSM. The mass of 125 GeV can be accommodated 
(with loop corrections giving 35 GeV) for certain range of values of the stop masses and left-right 
stop mixing parameter Xt = At — /itan/3. That range varies from Msusy = y/'^Wj^ ~ '''00 
GeV for the "maximal" stop mixing Xt ^ TeMgusY to MgusY ~ C(5TeV) for Xj = [3 E] 
Q Such values of the stop mass parameters are well consistent with the absence so far of any 
stop signal at the LHC but may look high compared to the standard expectations based on the 
naturalness arguments. Awaiting for more experimental progress, one may discard those, after 
all quite subjective expectations, or one may hope that a light stop is still hidden in the data, 
and investigate the ways of reconciling the 125 GeV Higgs mass with stop mass parameters 
below the values quoted above. This necessarily requires a beyond MSSM scheme, with a 
larger tree-level Higgs mass than in the MSSM. (It is useful to note that stop loop radiative 
corrections l^rn^^ = 25 (30) GeV can be reached with Msusy ~ 300 (400) GeV for the maximal 
mixing and with MgusY ~ 1-5 (3) TeV for Xt = 0). In this context, NMSSM has been discussed 
in the literature [9l |TOl [TTl [121 [IS]. In the pre-discovery era, NMSSM was discussed mainly as 
a scenario allowing for a Higgs mass significantly above the values predicted by the MSSM 
[m HSl Uni [TTI [TH] . The attention has been mostly focused on the new tree- level contribution 
to the Higgs mass coming from the singlet-doublet-doublet coupling in the superpotential, 
XSHuHd, which can be significant for low tan/3 values and 0{1) values of A. More recently, 
already after the discovery of the 125 GeV Higgs, the NMSSM has been discussed in the context 
of ameliorating the naturalness in the stop sector [T9l [20l [2T] also mainly for the same range 
of parameters. However, one may think that 125 GeV is close enough to the range expected in 
the MSSM, so that small corrections to the tree-level mass are worth considering. This is why 
it is interesting to investigate how significant may be the effect of the singlet-doublet mixing 
on the Higgs mass in the intermediate and large tan f3 region, where the MSSM tree-level value 
is ~ Mz- 

In NMSSM, there are three physical neutral CP-even Higgs fields, H^, H^, S which are 
the real parts of the excitations around the real vevs, Vu = wsin/3, Vd = f cos/3, Vg with 
■^^ = "^n + "^d ~ (174GeV)^, of the neutral components of doublets Hu, Hd and the singlet 
S (we use the same notation for the doublets and the singlet as for the real parts of their 
neutral components). It is more convenient for us to work in the basis {h,H,s), where h = 

^For a given value of the loop correction, the value of AfgusY depends also on other features of the SUSY 
spectrum, especially on the gluino mass and the stop mass splitting, the top mass and on unknown higher- 
order corrections to the Higgs mass. The value of AfsusY is particularly uncertain for AfsusY ^ 1 TeV, see 
e.g. Refs.iaiH]. 



HiicosP + HuSinP, H = H^sin/S — HuCos/S and s = S. The h field has exactly the same 
couplings to the gauge bosons and fermions as the SM Higgs field. The field H does not couple 
to the gauge bosons and its couplings to the up and down fermions are the SM Higgs ones 
rescaled by tan/3 and — cot/3, respectively. The mass eigenstates are denoted as s, h, H, with 
the understanding that h is the SM-like Higgs. 

In this paper we point out that s — h mixing effects can significantly contribute to ?ti/j, in- 
creasing it by 6-8 GeV, in the intermediate or large tan/3 region, making the NMSSM attractive 
also in that range of tan (5. Important for this scenario is the s — H mixing, which becomes 
significant for larger values of tan/3 and can lead to a suppression of the s ^ hh decay rate 
(where s is the singlet-dominated scalar), so that the LEP bounds on a lighter than 114 GeV 
scalar coupling to the Z boson based on that decay channel are evaded. |j The maximum mix- 
ing contribution to uih is then limited by much weaker bounds obtained from the s — > hadrons 
signature. The effects discussed in this paper require smallish values of the coupling A, C(O.l). 
Thus, the two regions of the NMSSM parameters, the low tan /3 one and the one considered 
here are clearly different. 

The s — H mixing at intermediate and large tan/3 has another interesting effect. It alters the 
decay rates of both s and h. They become correlated in an interesting way with the correction 
Am™'^ and show a variety of interesting patterns, depending on the lightest scalar mass. In 
particular, a striking signature of this mechanism can be a light scalar with strongly suppressed 
(enhanced) branching ratios to hh {gg, cc, 77) as compared to the SM Higgs with the same 
mass. The 77 decay channel is particularly promising for the search of such a scalar at the 
LHC. 

In section [2] we recall the structure of the CP-even scalar sector and discuss the effects of 
the s — h mixing on the Higgs mass, with the LEP bounds on a light scalar taken into account. 
In section |3J we discuss the potential role of the s — H mixing and the parameter range for 
which our mechanism can be relevant. In section |4] we give the predictions for the production 
and decays of the 125 GeV scalar if 5-8 GeV of its mass comes from the s — h mixing effects. 
In section |5] we briefly discuss the prospects for the discovery of a light scalar at the LHC and 
in section |6] we give a summary of the considered here scenario. 



2 CP-even scalar sector in NMSSM 

In this section we recall the necessary for us facts about the CP-even scalar sector of NMSSM 
[21]. Several versions of NMSSM has been proposed so far [251 [26l [27] . We would hke to keep our 
discussion as general as possible so we assume the NMSSM speciflc part of the superpotential 
tobeE] 

W^NMSSM = XSH^Hd + f{S) . (1) 



^The effects of the s — h doublet mixing in the CP-even scalar sector on the SM-like Higgs mass have also 
been discussed in the literature [52 US] for low tan /3 scenario. The mixing with the heavy doublet is usually 
ignored as it is small for small tan/3. 

^Explicit MSSM-like /i-term can also be present in the superpotential but it can always be set to zero by a 
constant shift of the real component of S. 



The first term is the source of the effective higgsino mass parameter, ^^s = ^Vg (we drop the 
subscript "eff" in the rest of the paper), while the second term parametrizes various versions of 
NMSSM. In the simplest version, known as the scale-invariant NMSSM, f{S) = kS^/3. 

We assume also quite general pattern of soft SUSY breaking terms (we follow the conventions 
used in [21]): 

-£soft 3 mljHu\^+mljHd\^+ml\S\^+{AxXHuHdS+-KA,S'^+mlH^Hd+-m'^S^+^sS+h.c.) . 

(2) 



Various versions of NMSSM studied in the literature [25], [261 HZ] belong to some subclass of the 
above setup. In the scale-invariant NMSSM, r/ig = m'g =^5 = 0. 
Let us parametrize the mass matrix of the hatted fields as follows: 



M' 



( Ml^ MIh Ml \ 

Ml^ Ml^ Ml, 
\ Ml Ml, Ml ) 



(3) 



where 



Ml- 
MIh 



M 



Ml cos^ (2/3) + {5mlY'''^ + \^v^ sin^ (2/3) 



{Ml 



A') sin" (2/3) + 



2Bii 



-At;2sin2/3 ( — 
2 \v. 



sin (2/3) ' 



M; 



hH 

2 
hs 



Mh 



(M|-AV)sin4/3^ 



M; 



Af(2/i-Asin2/3) 
= AwAcos2/3. 



(4) 
(5) 

(6) 

(7) 

(8) 
(9) 



Following [28], we introduced k = Ax + (9|/), while B = Ax + {dsf)/vs + ml/{Xvs) and 
T = {{dlff) + {dsfdlf) 



+ Af^KVg — —■ We neglected all the radiative corrections 



except those to Ml which we parametrize by {5ml)^'^'^. The first two terms in eq. (W| are the 
"MSSM" terms, with 



{5m 



2Nrad 
h) 



Sg'^mj 



iTC^TTl 



W L 



In 



ml ) M? 



SUSY 



X', 



12MgUgY 



(10) 



where Msusy = ^/'^ 



m.i,m 



t2 



{m^i. are the eigenvalues of the stop mass matrix at Msusy in the DR 
renormalization scheme) and Xj = Aj — /// tan /3 with At being SUSY breaking top trilinear 
coupling at Msusy- 

The third term in eq. (|4]) is the new tree-level contribution coming from the XSHuHd 
coupling. 

We recall that the eigenstates of M^ are denoted as s, h, H. We are interested in the 
parameter range such that m^ < m^h < m,H, so that the s — h mixing pushes the m^h up. We 
also require m^h < Im^g to avoid h -^ ss decays [29] . 



Quite generally, the mass of the SM-like Higgs reads: 

ml = M,\ + {6mlr'- . (11) 



The ((5?Ti^)™'^ term originates mainly from the ih,sj mixing and is positive (negative) when the 
singlet-dominated scalar is lighter (heavier) then the SM-like Higgs scalar. In the moderate and 
large tan/3 regime, the tree-level contribution coming from the XSH^Hd coupling is suppressed 
so one has to investigate in detail the potential effects of ((5?7i^)™^^. 

2.1 The effects of the s — h mixing on the Higgs mass 

In the case of no-mixing with if, the s — h mixing is determined by the 2x2 block of the mass 

matrix M^: 

'Ml, Ml\ 

ML Ml ; ' 

where the entries are given by eqs. (|4|), (pj) and (pj). The matrix M^ is diagonal in the basis 
s = gji + (3sS, h = v^l - gfh - v^l - (3^s. ^ 

In order to quantify the effect of the s — h mixing on the Higgs mass it is useful to introduce 
Amix such that: 

mh = Mhh + Amix • (13) 

Trading M|/j, Mf,, M^ for the two mass eigenvalues nih, rus and the coupling g^ of the singlet- 
dominated state to the Z boson (normalized to the corresponding couphng of the SM Higgs), 
one obtains a simple formula for A^ix: 

A^ix = m, - ^ml - ts {< - ml) ^ f (^'^ - ^) + <^(t) , (14) 

where in the last, approximate equality we used the expansion in ^^ <C 1. It is clear from 
the above formula that a substantial correction to the Higgs mass from the mixing is possible 
only for not too small couphngs of the singlet-like state to the Z boson and that m^ -C rrih 
is preferred. However, LEP has provided rather strong constraints on the states with masses 
below C(llO) GeV that couple to the Z boson because such states could be copiously produced 
in the process e"'"e~ — ;■ sZ. 

For those LEP searches that rely on the identifications of h and r in the final states |30j . 
constraints on g1 depend on the s branching ratios and the LEP experiments provide constraints 
on the quantity ^^ defined as: q] 

The LEP constraints on ^^^ are reproduced by the red line in the left panel of Figure ol 



^In the definition 15 it is implicitly assumed that the ratio r(s — > bb)/r{s —> tt) is the same as for the 
SM Higgs with the same mass. This is a very good assumption for NMSSM since the shh and stt couplings 
(normalized to the corresponding values of the SM Higgs) are the same at tree level. 



U.b- 












J 














'/ 












A. / 












f] 


1/ 










I 


/ 










A A 


K/^' 


^ 




U.Z" 




A/^ 


Ayv 


J 






rt 1. 


. r 


/ 


/ 














J 








0- 


■-^ 


' 


' 


1 


— ' — 





60 



.max 
A . 5 
mix 



70 



80 



90 



100 



110 




m 



70 80 90 100 110 

m 



Figure 1: Left: The 95 % CL upper bounds on ^^^ (the red hne) or $,jj (the green hne). The red 
hne was obtained using the observed hmits presented in the column (a) of Table 14 in Ref. |30j . 
while the green line corresponds to Figure 2 of [31]. Right: The LEP limits translated to the 

assuming ^j = g^ and 6 



upper limits on Amix using eq 
respectively. 



14 



■'jj 



gj. for the red and green line, 



Since we assume in this subsection that s mixes only with h, all the couplings of s are those 
of the SM Higgs multiplied by a common factor g^. This implies that the branching ratios of s 
are exactly the same as for the SM Higgs. Therefore, the limits on ^^^ depicted by the red line 
in Figure ^ are, in fact, also the limits on g'^. Using eq. (14) we can translate the constraints 
on gl into limits for the maximal allowed correction from the mixing, A™?*^, as a function of 
rris- These are presented in the right panel of Figure [Tj Notice that in this case the correction 
from the singlet-doublet mixing can reach about 6 GeV in a few-GeV interval for rus around 95 
GeV, where the LEP experiments observed the 2a excess in the bb channel. This is interesting 
since such correction combined with the tree-level values ~ Mz (for moderate and large tan/3) 
gives rrih ~ 125 GeV with Am]^'^ ^ 30 GeV. However, for m^ < 90 GeV the allowed value of 
A™fjf drops down very rapidly to very small values. 

In Figure [2J we present an example of the NMSSM parameters for which mixing with H is 
negligible and Amix ~ 6 GeV can be obtained. Note that the h — s mixing, thus also Amix, 
grows with tan /3 as a consequence of the suppression of the second term in the parenthesis in 
Mlg at large tan/3, see eq. (|8). This example demonstrates also the fact that A is generically at 
most C(O.l). Larger values of A typically lead to too large M|^ (after taking into account the 
LEP limit on the chargino mass which imply /z > 100 GeV) leading to a negative determinant 
of the mass matrix. Therefore, this scenario is the most natural at moderate and large tan/3, [j 



^In the scenario with Amix > 0, values of A '^ 0.6 that lead to substantial tree-level contribution to rrih at 
small tan/3 can only be obtained if (2/i — Asin(2/3)) (which enters Af^^) is finely-tuned to be below ©(lOGeV). 
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Figure 2: Amix [GeV] (black solid line), g\ x 10 (green solid line) and ^ x 10 (red solid line) 
as a function of tan /3. The remaining parameters are fixed to be: vnh = 125 GeV and /i = 150 
GeV, m, = 95 GeV, thh = 1500 GeV, A = 1200 GeV and A = 0.08. The constraint on ^^ x 10 
is depicted by the dotted red line. The regions with the solid red line below the dotted red line 
are allowed by the LEP data. 



It is also clear from Fig. IT] that similar correction 0{5) GeV to the Higgs mass can be 
obtained from the s — h mixing for a larger range of the singlet-dominated scalar mass tti^, 
provided one can evade the LEP bounds given by the red curve in the left panel of Fig. [T] by 
suppressing the shb and srf couplings. This is because in such a case, s decays predominantly 
into charm quarks and gluons and 6-tagging cannot be used to enhance the signal over back- 
ground ratio so the most stringent constraints on IJl come from the flavour independent Higgs 
searches in hadronic final states at LEP [Slj. Those searches give constraints on a quantity si- 
delined as: 



e^, =^^xBR(s^jj) 



(16) 



which are reproduced by the green line in the left panel of Fig. [T} Noting that for suppressed 
sbb and srf couplings, BR(s — )■ jj) ~ 1 so ^^ 



g1, we can translate those constraints into the 
upper bound on Amix. Indeed, the upper bound A™fx is then given by the green curve in the 
right panel of Fig. [Tj 

We show in the next section that s — H mixing can significantly change the decay rates of 
s and also of h. H Firstly, the A™?f shown by the green line in Fig. Q^ can then be obtained 
for a broad range 60 GeV < m^ < 110 GeV, and secondly the decay rates of s and h can have 
interesting patterns. 



^ A suppression of the s ^ bb decay rate is possible for any value of iris but for nis in the few-GeV interval 
around 95 GeV there is no gain in A^f^ because the red and green curves in Fig. fll practically overlap there. 
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3 Singlet mixing with both doublets and the suppression 
of the sbb coupling 

We now go back to the general case in which mixing with H may be present. As in the previous 
section, we begin with the imphcations following from the general structure of the mass matrix. 
Mixing with H leads to the modification of Higgs couplings to fermions. Denoting the mass- 
eigenstates s, h, H hy x = g^h + pi H + l3x s we get 

a. =^. + /3fHan/3, (17) 

a=^.-/3f)cot/3, (18) 

Cv. = 9. , (19) 

where x is s, h or H. Note that the couplings to the vector bosons depend only on the h 
components, as in the case of only {h, s) mixing discussed in the previous subsection. 

In the region of moderate and large tan/3 even small component of H in the singlet- 
dominated Higgs may give a large contribution to the couphngs to b quark due to tan/3 en- 
hancement. On the other hand, the couplings to the up-type quarks are almost the same as 
those to the gauge bosons, Ct^ ~ Cy^- Particularly interesting is the case when g^ has the 
opposite sign to f3s because then Cb^ <^ Ct^ , Cy, is possible. In the regime Cb^ <^ Ct^ , Cy^ , the 
(otherwise dominating) s branching ratios to bb and rf are strongly suppressed and s decays 
mainly to gg and cc. The ratio r(s — )■ gg)/T{s — )■ cc) is roughly the same as for the SM Higgs 
so e.g. for TUs = 90 GeV it equals about 1.5 [321 [33] and approximately scales like m"^ for other 
masses [M]. In this regime the standard LEP Higgs searches [30] that used 6-tagging cannot be 
applied to constrain this scenario. In such a case, the most stringent constraints comes from the 
flavour independent search for a Higgs decaying into two jets at LEP [31j. These constraints 
are weaker and allow for values of g^ above 0.3 for mg around 100 GeV and the limit rather 
slowly improves as rris goes down, as seen from the left panel of Figure [Tj In consequence, 
the constraints on A™|^ are also weaker. As can be seen from the right panel of Figure [Tl 
when s -^ bb decays are suppressed Amix above 5 GeV is viable for a large range of iris with a 
maximum of about 8 GeV for nig around 100 GeV. 

We should also comment on the fact that for C^^ <^ Ct^ , Cy^ the s branching ratios to the 
gauge bosons are also enhanced (with respect to the SM Higgs predictions) by a factor that can 
exceed 10. In spite of such a large enhancement the Higgs searches in these channels performed 
at LEP [35] are less constraining than the above-discussed searches with hadronic decays. On 
the other hand, the LHC searches in the diphoton channel may, in principle, have a potential 
to give additional constraints on this scenario (i.e. reduce the allowed value of A™?f). In fact, 
s — )■ 77 decays could already be seen at the LHC but the SM Higgs searches in the diphoton 
channel have been performed only for masses above 110 GeV. This will be discussed in more 
detail in section |5] 

In the above discussion we assumed that the strong suppression of the s coupling to b quarks 
is possible. Let us now discuss in which part of the NMSSM parameter space such situation 
may hold. As already stated, this may happen only for not too small values of tan/3 and a 
negative ratio (3s /~gs- This ratio can be expressed in terms of the M^ entries and vfis in the 



following way: 



^ = —^. i :—■ (20) 



At large tan/3, M^^ ~ — 2(M| — A^f^)/tan/3 is very small so the second terms in the numerator 
and the denominator are typically sub dominant H which means that /3s /gg is negative if 
Mjj^Mlg < (we recall that m^ < M|^ in our case) which leads to the following condition for 
the NMSSM parameters: 

A(/itan/3-A) >0, (21) 

which is satisfied only if /lA > 0. In the following discussion we will assume, without loss of 
generality, A > and /i > 0. It is straightforward to show in the limit of large tan/3 that Cb^ 
may vanish only if 



2 



where 



9 2A 
^ >-r^ (22) 



r^ ^ Mh^< . (23) 



If the condition (|22|) is satisfied then Ci,^ ~ corresponds to two values of tan /3: 

(24) 




which in the limit r^ ^ ^ are given by: 



A / A2 \ ^ 2/ir2 

- 1 + -^^ V tan/3 ^^ 



tan/3^-(l + — -^) V tan/3^^. (25) 



Let us now demonstrate some numerical examples in which the suppression of Cb^ is present 
and substantial values of A^m is obtained without violating the LEP constraints. In Figure [3] a 
tan /3-dependence of A^ix is presented. It is clear from this Figure that substantially larger Amjx 
is consistent with the LEP data due to the suppression of the sbb coupling. In the left panel, 
rris = 100 GeV and Amix can be almost 8 GeV. The role of the suppression of the sbb coupling 
is even more important for lighter singlet-dominated states. In the right panel, m^ = 75 GeV 
and Ajnix can reach 6 GeV at large tan/3, while without the suppression it would be below 2 
GeV. 

It can also be seen from Figure [s] that there exist values of tan (3 for which the sbb coupling 
is strictly zero. Nevertheless, such a strong suppression is not necessary to avoid the LEP 
constraints. In fact it is enough to suppress BR(s — )• bb) by about 25% in the case of rus = 100 
GeV and by a factor of three for mg = 75 GeV. This implies rather large range of tan (3 with 
significant correction from the mixing consistent with the LEP data. 

'^Strictly speaking, the second terms in the numerator and the denominator can dominate for A — > because 
then Mfjg -^ 0. In such a case (31 /g^ is negative if M'^^j > which is possible only if \^v^ > M|. However, 
for X^v^ > M"^ and A — j' the mass matrix has a negative eigenvalue if M^^, < Af,-^^ (which is a necessary 
condition for s to be lighter than h). 
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Figure 3: Amix [GeV] (black solid line), g\ x 10 (green solid line) and ^ x 10 (red solid line) 
as a function of tan/3. For easy reference the LEP constraints are depicted by the dotted lines 
with the same color coding as in Figure [1} A point in parameter space is consistent with the 
LEP data if the red solid line is below the red dotted line and the green solid line is below 
the green dotted line. In both panels vn}^ = 125 GeV and fi = 150 GeV. In the left panel: 
rris = 100 GeV, m^ = 500 GeV, A = 600 GeV and A = 0.06. In the right panel: m, = 75 GeV, 
ruH = 1000 GeV, A = 800 GeV and A = 0.08. 



In our analysis we use the eigenvalues of the Higgs mass matrix as input parameters while 
the diagonal entries of this matrix are output parameters. Such procedure is justified because 
any values of the diagonal entries can be obtained by adjusting the soft terms in appropriate 
way. However, it is natural to ask whether the required values of soft terms are reasonable. One 
cannot answer this question in a model-independent way so let us focus on the no-scale version 
of NMSSM which is the most popular one and calculate the soft terms in some representative 
examples. In such a case, /i = Xvg, B = Ax + nvg and A = Ax + 2/tt>s. Requiring the correct 
electroweak minimum, for the parameters used in the left panel of Figure [3] one obtains for 
tan/3 = 25 (corresponding to Amix ~ 7.3 GeV): 



A^ = -2111 GeV, Ax 



-467 GeV, 



K 



0.213, 



m 



Hu 



-(162 GeV) 2, m^ = (480 GeV)^ 



TTlc 



(109 GeVf 



(26) 



while for the parameters used in the right panel of Figure M one obtains for tan /3 = 35 (corre- 
sponding to Ajnix ~ 6 GeV): 



A^ 



-2427 GeV, A, 



-419 GeV, 



K 



0.325, 



TJl 



\^ = -(161 GeV)^ ml^ = (990 GeV)^ 



rric 



(85 GeV)^ 



(27) 



From the above two examples it should be clear that large values of Amix can be obtained for 
rather natural values of the soft parameters. In particular, the values of k can be consistent 



10 



with the upper bound, Kmax ~ 0.65 (for A < 0.1) [21], from the requirement of perturbativity 
up to the GUT scale. 



4 Production and decays of the 125 GeV Higgs 

The mixing effects affect not only the branching ratios and production cross-section of s but 
also those of h. Moreover, they are correlated so the scenario may be tested also by the 
measurements of the signal strengths for the 125 GeV Higgs. In order to set a notation let us 
define the signal strengths modifiers as: 



,(h) ^ O-JPP ^ h) X BR{h -^ i) 

a^^ipp ^h)x BR^^{h -^ i) 






In the case of the h — s mixing, with the effects of H neglected, all the h couplings are 
multiplied by a common factor a/1 — 'gl- This implies that all the h branching ratios are the 
same as for the SM Higgs while the production cross-section (in all channels) is smaller by a 
factor 1 —~g1 so Ri = 1 — ~g1 for all channels. This means that, after taking into account the 
LEP constraints, A^i^ > 5 GeV implies 0.75 < /?f ^ < 0.83. 

In the full 3x3 mixing case at large tan/3, the couplings to the up-type quarks are almost 
the same as those to the gauge bosons, Ct,^ ~ Cy^ = V'l — gl so the production cross-section is 
still smaller than the SM prediction by a factor l—g^- However, the couplings to the down-type 



fermions can be substantially modified, as seen from eq. (17). Since the bb channel dominates 
the decays of the 125 GeV SM Higgs (BR(/iSm ^ 5^) ~ 53% ^nd BR{h^^ -^ rf) ^ 6%) such 
modifications lead to important effects for all the other branching ratios. If f3j^ /gt i^ negative 
(positive) then the h couplings to b and r are smaller (larger) than in the SM which leads to the 
enhancement (suppression) of the Higgs branching ratios to the gauge bosons and two photons. 
This ratio is given by 



/3f ) _ _ Mf,, (ml - Ml) + MlMs 
~^' {mIh - mi) [ml - Ml) + (m^. 



(29) 



and its sign is: 



sgn I ^1 = -sgn I Mi, + , ^'^^'^^ , I . (30) 



9h 



ml 




Since Mfj^ is small, the enhancement (suppression) of the h coupling to b requires Mfj^M^^ < 
(> 0). Note that this is the opposite condition to that for the s coupling so if the s coupling 
to b is enhanced (suppressed) then the h coupling to b is suppressed (enhanced). As it was 
discussed in the previous section, for the m^ in the range between about 90 and 105 GeV, 
Amix can exceed 5 GeV with the LEP constraints satisfied independently of the sbb coupling 
and both discussed above options are interesting. From the current experimental viewpoint the 
suppressed h coupling to b is more welcome in order to compensate the suppression of the h 
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production cross-section and end up with Ryy ~ 1 (where ^ = ly or Z). H However, given the 
present tension between the CMS and ATLAS results the case with the enhanced h couphng 
to b is certainly not excluded. 

The predictions for R\^' are very similar to Ryy because the reduced couplings to top and 
W (which contribute to the /i — )■ 77 decay in the SM) are almost the same at large tan/3, 
Ct^ ~ Cv^. The enhancement of -R77 over Ryy, which is preferred by the ATLAS data, is 
possible only if contributions of SUSY particles to the /i — )■ 77 decay width is non- negligible. It 
was shown in [28] that such enhancement can be substantial for light higgsinos and A ~ 0{1). 
However, we found that this effect is small in our case since A is required to be (9(0.1) at most. 
The most promising way to obtain the 77 enhancement would be the presence of very light 
staus with strong left-right mixing which may be possible if tan /3 is large |37] . 

4.1 s with strongly suppressed couplings to b and r 

It is crucial to note that the couplings of h and s to 6 are correlated. It is the purpose of this 
subsection to investigate the implications of the strongly suppressed s couplings to b for the 
production rates of h. 

In order to study quantitatively the correlation between the correction to the Higgs mass 
from mixing and the production rates for h we performed a numerical scan over the NMSSM 
parameter space for various values of rris and uih while keeping fixed rrih = 125 GeV. In the 
scan we also fixed fi = 150 GeV. For other values of /i the results of the scan are the same 
provided that the following transformation of parameters is used: 

fx-)- kn, \-^ \/k , A -^ A;A . (31) 

This is because M^^ and Mfjg are invariant under the above transformation while M^h i^ 
only marginally affected so its impact on the numerical results is negligible. The remaining 
parameters where scanned on a grid, see Table [T] for the scanned parameters ranges and step 
sizes. In order to emphasize that obtaining substantial values of Amix does not require any 
fine-tuning the grid is not dense, as clearly seen from Table [Tj In Figure |4] a scatter plot 
of Amix versus rris is presented. The LEP constraints discussed before have been taken into 
account. It can be seen that Amix up to about 9 GeV can be obtained for m^ ~ 100 GeV but 
for such large values of A^ix Ryy < 0.5 is predicted, which is in tension with the LHC Higgs 
datajj Nevertheless, demanding Ryy > 0.5, Amix about 8 GeV can be reached. Notice also 
that Amix ^ 5 GeV with Ryy > 0.7, which is well consistent with the LHC data within the 
experimental errors, can be obtained for a wide range of values between m/i/2 and 105 GeV. 

The reduction of Ryy is due to the h production cross-section suppressed by a factor 1 — g^ 
and the suppressed BR(/i — )■ VV), as a consequence of the enhanced hbb coupling. However, 
for rus between about 90 and 105 GeV, where the LEP constraints on g^ are not so strong 



^Such a scenario was investigated for low tan /3 in |36| with a special attention to possible 77 rate enhancement 
for the 125 GeV Higgs. _ 

I are slightly larger than the corresponding 
BR(s -^ jj) = 1, i.e. gl = ^^-j is assumed, 



^Notice that maximal values of Amix ior a given rrig in Figur 
values in the right panel of Figure 1 This is because in Figure 



while for the points from the numerical scan that are consistent with the LEP data the sbb and stt couplings 
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are not exactly zero so BR(s — )• tt) > leading to g^ > ^^ 



■g 4 : 



3 
2 
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Figure 4: Results of the numerical scan presented in the A^ 

correspond to different values of Ryy. 

for the red, green, black and yellow points Ryy is larger than 0.5, 0.7, 0.8 and 1, respectively. 

The points with larger values of Ryy are overlaid on the points with smaller Ryy. All the 

points satisfy the LEP constraints. 



-nis plane. Different colours 
The blue points are characterised by Ryy < 0.5 while 





rriH [GeV] 


A 


A [GeV] 


tan/3 


Minimal value 


250 


0.05 


100 


10 


Maximal value 


2000 


0.15 


3000 


60 


Step size 


250 


0.01 


100 


5 



Table 1: The parameter ranges and step sizes used in the numerical scan. 



and rUg is still significantly below 125 GeV, suppression of the sbb coupling is not necessary for 
obtaining substantial values of Amix. Therefore, in that range Ryy > 1 can be obtained with 
Amix ^ 5 GeV. Such solutions are characterised by the enhanced sbb coupling and suppressed 
hbb coupling. 



Since for moderate and large tan/3, O 



{h) 

■77 



Ct, the predictions for R. 



is almost the same 



ouiue loi uiuueiaie anu laige lan/j, L/y ><:- c/^, iiie pieuicLions lui rt^-y ib aiiiiobi ine same 
as for Ryy, presented in Figure 4l After taking into account the higgsino contribution to the 



becomes slightly larger than Ryy. 

Ah) 



For /i = 150 GeV, which was used 



/i — )• 77 decay rate R. 

in the numerical scan, R)y!y' is typically enhanced with respect to i?^'^ by a few percent. 

We should stress that our analysis is performed at tree level. It is well known that at large 
tan/3 SUSY threshold correction to the bottom quark Yukawa coupling may be substantial 
(38| |39] . If those corrections act in such a way that the loop-corrected C^^ is smaller than the 
tree-level value then the h branching ratio into gauge bosons is enhanced. Thus, in principle 
some regions of the NMSSM parameter space may exist in which Amix reaches 8 GeV and Ryy 
is around one, without violation of the LEP constraints. However, a detailed study of such 
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corrections is beyond the scope of this paper. 

5 Prospects for discovery of s at the LHC 

Let us now discuss prospects for discovery of s at the LHC. What the experiments observe is 
the product of the production cross-section and the branching ratios: 



j(s) ^ ff(pp ^ s)x BR{s -^ i) 

aSM^pp ^h)x BR^^{h -^ i) 






If the s branching ratio to bb is not strongly modified as compared to that of the SM Higgs, the 
signal strengths in all channels are universally suppressed i?j ~ g^. For nis in the range 90-105 
GeV, where Amix > 5 GeV is possible without strong sbb coupling suppression, R\^ < 0.25. 



^ In that range of rris the LHC experiments have the best sensitivity in the s -^ bb decay 
channel (the 77 and tt channels may also be relevant, especially for rris > 100 GeV). The LHC 
experiments do not provide limits, nor the expected sensitivities, for the masses below 110 GeV 
in their searches for the SM Higgs. However, the sensitivity of the search in the bb channel very 

(s) 

weakly depends on nis in this range so one can estimate that the expected sensitivity to -R^g 
is about 0.9 with the data that have been analysed so far i.e. 5 fb~^ of the 7 TeV data and 13 
fb~^ of the 8 TeV data [HI |12]. From a naive extrapolation to higher luminosities one expects 
that about 200 fb~^ of the 14 TeV run will be required to test this scenario. 

5.1 s with strongly suppressed couplings to b and r 

It should be clear from the previous section that the scenario with a strong suppression of 
the sbb and stt couplings can be constrained by the precision measurements of the 125 GeV 
Higgs couplings. Even more interesting is the fact that the LHC is already well prepared for a 
discovery of s. This is because in this scenario the total decay width of s is strongly reduced so 
all the s branching ratios, except those for the s decays to the down-type fermions, are strongly 
enhanced. 

Particularly interesting is the 77 final state. [^ In Figure [s] we present the predictions for 
i?77 assuming maximal value of g'j. consistent with the LEP s — > jj data (corresponding to 
maximal value of Amix allowed by the LEP data) as a function of rus- In the extreme case when 
the sbb and stt couplings are suppressed to zero (the black line in Figure [5|, the 77 signal from 
s decays is stronger than that from the SM Higgs with the same mass for the whole range of 
rris- For rris around 100 GeV the enhancement can almost reach a factor of three. 

As we already mentioned, it is not necessary to suppress the sbb and stt couplings exactly 
to zero. In fact, it is enough to suppress them to the level for which the LEP constraints on 



^fr are satisfied. The blue line in Figure 5 correspond to the minimal suppression of sbb and 
STT couplings required to satisfy the constraints on ^^^. Even in this case R^^j > 1 for a wide 



^"Particularly interesting possibility is the singlet-like Higgs with mass about 98 GeV because it can explain 
the LEP excess in the bb channel [40] , 

^^The possibility of large 77 rate enhancement for the singlet-like NMSSM boson was noticed in Ref. ^^. In 
contrast to the present paper, in Ref. [43] small values of tan/3 sa 3 were considered. 
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Figure 5: The predictions for Wy^.j assuming maximal value of g^ consistent with the LEP 
s — )■ jj data as a function of rris. The black line corresponds to the case when the sbb and stt 
couplings are suppressed to zero. The blue line correspond to the minimal suppression of the 
sbb and stt couplings required to satisfy the LEP constraints on C,?t. 



range of m^ between about 60 and 90 GeV, and around 100 GeV. Small values of Wf'-y are 
possible only in the few-GeV interval around 95 GeV where the LEP limits on C,?t and ^? are 



66 



^jj 



comparable. 

Despite such significant enhancement of the 77 rate the LHC experiments do not constrain 
this scenario for rris below 110 GeV because the Higgs data have not been analysed in that 
region. For rris = 110 GeV, the current observed CMS upper limit [H] (based on 5 and 20 fb~^ 
of the LHC data at 7 and 8 TeV, respectively) on R>f^.j is about 0.6 which already constrain the 
allowed values of g^, thus also A™|^, for this particular mass. Therefore, one can expect that 
the LHC searches are sensitive enough to probe this scenario in the 77 channel also for smaller 
values of m^. 

Since the expected limit on Ky^ for rus = 110 GeV with the current data is about 0.6 |44j . 
and the sensitivity gets worse quite slowly when the mass goes down, a naive extrapolation of 
the available analyses suggests that the LHC could have already set the limits on Wf^-j ~ (9(1) 
for masses below 100 GeV using the available data if these were analysed. 



6 Conclusions 



We have studied in detail the mixing between the three physical scalars s, h and H of the CP- 
even scalar sector of the NMSSM. In a large parameter range, it can lead to several interesting, 
often correlated, effects. First of all, the s — h mixing can give 6-8 GeV contribution to the 
mass of the SM-like scalar h in the moderate and large tan/3 region and with A ~ (9(0.1). This 
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is interesting because the 125 GeV mass is then obtained with significantly lower stop masses 
in the stop-top loops. The geometric mean of the stop masses, Msusy, can be below about 
400 GeV (2 TeV) for the maximal contribution from stop mixing (with no stop mixing at all). 
Thus, the NMSSM is interesting also beyond the usually considered region of low tan /3 and 
A~C>(1). 

The s — h mixing contribution to rrih depends mainly on the mixing angle between the two 
fields (i.e. on the sZZ coupling g^) and on their mass difference. Thus the effect is constrained 
by the LEP bounds on g^ vs rris obtained assuming for s the SM Higgs branching ratios and 
with b and r identification in the final state, or without such particle identification assuming 
BR(s — )■ hadrons) = 1. The two experimental bounds almost overlap for rris in the 5 GeV- 
interval around 95 GeV but for other s masses the bound based on the detection of two non- 
identified hadronic jets is much weaker. For rris in the 5 GeV-interval around 95 GeV the 
(9(6 GeV) mixing contribution to rrih can thus be obtained independently of the decay modes 
of s. However, for other values of rrig the mixing contribution is much smaller, unless the LEP 
bound based on the b and r identification is evaded, i.e. if s — )■ 66 is suppressed strongly enough. 
In the latter case, the 5-8 GeV effect is obtained for the range 60-110 GeV of rris, consistently 
with the LEP bound based on the search for two hadronic jets, with BR(s — )■ hadrons) = 1. 

Interestingly enough, a strong s — )> 66 suppression can be present due to the s — H mixing 
(with negligible effect on the s — h sector), which is important in the considered region because 
of the tan /3 enhancement of the scalar down quark couplings. Thus the LEP bounds can be 
evaded. The lightest scalar s has then enhanced branching ratios into ZZ*, WW* and 77. The 
latter one is a particularly promising signature for the LHC searches for a scalar lighter than 
110 GeV, with suppressed 66 decay channel. The signal strength in the 77 channel of this scalar 
may be larger than that of the SM Higgs, even by a factor of three. In fact, if such singlet-like 
scalar with mass below 110 GeV really exists it could have already been discovered at the LHC 
if the already collected data were analysed in this range of masses. Thus, we strongly encourage 
the ATLAS and CMS collaborations to extend their Higgs searches in the 77 channel to masses 
in the 60-110 GeV range. 

The s — H mixing modifies also the h decays, in a way anti-correlated with the s decays. The 
ones suppressed for s are enhanced for h and vice versa. Thus, the large mixing contribution 
to rrih can be present together with a variety of interesting patterns for the h production and 
decays. If rris is between about 90 and 105 GeV, the mixing correction to rrih exceeding 5 
GeV does not require the suppression of the s66 coupling and e.g. BR(/i — )■ 77) can be either 
enhanced or suppressed as compared to the SM prediction. If rris is smaller or larger than the 
values given above, the large mixing effect is generically correlated with suppressed rates in 
ZZ, WW and 77 channels and enhanced ones in 66, tt channels for h. The magnitude of that 
suppression (enhancement) depends on the particular choice of parameters. 

The effects considered in this paper do not require any particular fine tuning of the NMSSM 
parameters and are present in a large part of parameter space. 
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